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Use the XOR function to implement modulo 2 division 
when calculating cyclic redundancy checks 



Recently I needed to im- 
plement the XMODEM 
cyclic redundancy check 
(CRC) option on my 
XCom9 modem pro- 
gram. Despite the many 
programs available for 
calculating CRCs, I had some difficulty 
understanding the math behind the calcula- 
tions. The details became clear after much 
research and experimentation. 

The starting point for all CRCs is fancy 
linear algebra. The CRC is defined in terms 
of message polynomials, generator poly- 
nomials, and so on. As Perez. Wismer, and 
Becker state in "Byte-wise CRC Calculations" 
IIEEE Micro, lune 1983). 

In a system employing CRC the message be- 
ing transmitted is considered to be a binary 
polynomial MIX). It is first multiplied by X k 
and then divided imodulo 2 1 by an arbitrary 
generator polynomial G/X) of degree k 
which results in a quotient OIXI and a 
remainder. . . 

It sounds confusing, but if you can under- 
stand how to apply the math, you can 
design more efficient programs and spot 
many erroneous ones. 

THE "MATH-NESS" TO THE METHOD 

The design of the polynomial CIX) is ex- 
tremely complex. You need to pick one that 



produces CRCs that are good at detecting 
errors. Fortunately, many G(X)s exist 
already. Table 1 contains the two most com- 
mon CIX)s in an 8-bit-byte environment. 

Let's calculate the CRC for the letter T. 
0101 0100 in binary. MIX) is the message 
as it is transmitted. You transmit a charac- 
ter's least significant bit (LSBI first, so MIXI 
becomes 00101010. Then you divide 
modulo 2 as shown in figure 1. (Modulo 2 
means you use the XOR instruction instead 
of the normal add and subtract.! Work it 
through according to the process shown in 
figure 2. Note that the CRC result is given 
in reverse order, that is. most significant bit 
IMSB1 on the right. LSB on the left. 

Bit-oriented Algorithms 

Using the long-division approach, if you had 
only a single zero bit to send, you would 
get the result shown in figure 3. If you had 
two bits to send, a zero and then a one. long 
division would produce the result shown in 
figure 4. The first remainder in figure 4 is 
the same as the first remainder in figure 3 
except that its LSB has been XORed with 
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Tahle I' Tfif am most common CiXIs [generator polynomials) in an §*t-%te 
J™ ronment Y ol can code the CtX, as a 17-bit binary, or hmdecmal number. 
TheZmdZos represent the coefficients of the different powers of X m * 
polynomial. 

CRC-16 (Bisynchronous) 

U x „ + x> + 1 = 1 1000 0000 0000 0101 

CRC-CCITT (SDLC, X.25. XMODEM) 

X" + X" + X s + 1 = 1 0001 0000 0010 0001 
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Fieure 1- The initial modulo 2 dmsion used to calculate the CRC for the letter T. 
Hote that the letter is given LSB first; in other words, it isjeversed. 
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Fieure 2 The entire long-division calculation process used to derive the CRC for 
2 Tetier T The apparent subtractions in this process are XORs: therefore there are 
no %SmC A ' represents a urn that has no further part in the calculation, 
and Represents a zero that was "brought down from the dividend I,* 
Therefore the line 00 1000 0101 0010 100 represents a remainder of 00 000 
0010 10 Note that the CRC is aiven with the LSB on the left. Thus the 
CRC for the L T = 54 is (MSB jKl 000, 0.00 .0.0 000.. or I4A1 
hexadecimal. 



the second data bit. This similarity will 
always be true because of the way the 
MIX) polynomial is built. 

This observation leads to the follow- 
ing bit-by-bit algorithm for calculating 
the CRC: 

1. Write down the first data bit (zero 
or one) to be transmitted. 

2. Write down 16 zeros to its right. 

3. Divide the 17-bit number by the 
17-bit CRC polynomial using 
XOR instead of subtraction. 
Make a note of the remainder, 
which is the CRC. 

4 Get the next data bit. 

5 XOR this bit with the LSB (left- 
most bit) of the CRC in step 3. 

6. Append a zero to the right-hand 

end of the result in step 5. 
7 Divide the 17-bit number from 
step 6 by the 17-bit CRC poly- 
nomial. Use XOR instead of sub- 
traction. The remainder is the 
CRC. 

8. Repeat steps 4 through 7 until 
there are no more data bits. 
(You can replace steps 1 through 
3 with a single step to initialize 
the CRC to zeros.) 

Thus, you can calculate the CRC a 
bit at a time, for any number of bits, 
with almost no extra calculation over- 
head. In summary, the steps involved 
in the long-division method are 

1 The message bits are written 
down in the order in which they 
are transmitted, from left to 
right, that is. LSB on the left. 

2 Sixteen zeros are appended to 
the right-hand end of the binary 
number formed in step 1. 

3 The generating polynomial is 
written down MSB first, that is. 
on the left. 

4. The division is done modulo 2, 
that is, using XOR instead of 
normal subtraction. 

5. The CRC is the remainder after 
all data bits have been pro- 
cessed. The LSB is on the left. 

Figure 5 translates these steps into a 
flowchart. 

HARDWARE IMPLEMENTATIONS 

One disadvantage of the flowchart in 
figure 5 is that it requires a 17-bit reg- 

(HMItifluM 
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Figure 3: The long-division CRC calculation for a single mo bit. 
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Figure 4: The long-division CRC calculation for the two bits 01. 
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Figure 5: A flowchart representing the steps involved in the long-division method of 
calculating a CRC. 
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A powerful computation and 
documentation tool for your 
IBM-PC. 

With MathCAD you simply and inter- 
actively create, edit and display for- 
mulas on the screen the way you are 
used to writing them. Equations are 
instantly computed and the results dis- 
played on the screen as a single num- 
ber or a plot. Text may be added to the 
screen and everything may be printed 
out as an integrated document. 
MathCAD has built-in hyperbolic and 
circular functions, performs all calcu- 
lations with real and complex num- 
bers, performs iterative calculations, 
handles all units, performs error 
checking and dimensional analysis 
and much more... 

The price of MathCAD-$189. 

In Massachusetts add 5% sales tax. 

To order send check, p.o., call us 
with your MasterCard number or 
call us for the nearest dealer. 

I.SOO.MathCAD or 617.577.1017 

Math Soft 



One Kendall Square 
Cambridge. Massachusetts 02139 
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Figure 6: T(ie division process ij Q - 1 
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Figure 7: The division process if Q = 0. 
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Figure 8: The division process regardless o\ the value of Q. 
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Figure 9: The process shown in figure 8 translated into a hardware circuit using shift registers and XOR gates- [The boxes 
are stages to a shift register: the shift register is only 16 stages long ] Note that the XOR gate going into R15 is superfluous, 
since Q XOR is always Q. Note also that in the long-division method, the calculations are done before "bringing down" the 
next zero bit. Similarlu here, the XORs are performed before the shift. 
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Figure 10: The reverse drawing of figure 9: the LSB is on the right instead of on the left. 



Ill BYTE * SEPTEMBER 1 986 




CALCULATING CRCS 



ister. However, after the division is 
complete, the left bit will always be a 
zero. Either it was a zero to begin 
with, or if it was a one. it was XORed 
with the X'« bit of the polynomial 
which is also a one, producing a zero 
Thus, you can XOR with only the 16 
LSBs of the polynomial. In the case 
of the CRC-CCITT polynomial. X" + 
X" + X'+ 1 (I OOOI 0000 0010 0001 
m binary or 1 1021 in hexadecimal) 
you can XOR with 1021. 

Let the bits of the CRC register be 
RO. R1 . up to R15: let the data bits be 
DO up to D7; let the polynomial bits 
be XO up to X15; and define Q = D 
XOR RO. that is. set Q equal to the 
current data bit. D. XORed with RO 
In every case, the LSB is bit 0. 

Then, if Q - I. the division process 
looks like that in figure 6. If Q = D 
XOR RO is zero, then the division 
reduces to that in figure 7. But when 
Q - I. the new R15 is always I and 
you XOR R4 and R11 with X12 and 
X5. which are 1. If Q - o. then the 
new R15 is zero and you XOR R4 and 
R11 (and all other Rns) with zero. 

You can combine the two cases as 
shown in figure 8. Hirning that pro- 
cess into a circuit using shift registers 
and XOR gates (see figure 9) is 
straightforward. You can also draw the 
circuit so that the LSB (RO) is on the 
right (see figure 10). Both diagrams 
are convenient depending on the kind 
of software algorithm to be derived. 
Note that in hardware the XOR is 
done before the shift because of the 
way flip-flops work. The important 
point is that the new R10 - (old R11) 
XOR (old RO) XOR (new data bit). 

Bit-by-Bit Software 
Algorithms 

When you do the calculations in soft- 
ware, you don't have to do the XOR 
before the shift. The important thing 
in software is to avoid having to deal 
with 1 7-bit values that don't fit into a 
variable. A software routine also 
doesn't need to use the Q = D XOR 
RO result to drive a gate. You can pro- 
gram the polynomial directly into the 
code as a constant. 

You can now derive a software rou- 
tine based on figure 10 (or figure 9. 
in which the order of storing the bits 
of the CRC is reversed). If you assume 



that the CRC result is kept in a 1 6-bit 
integer with R15 (the MSB of the CRC) 
in the high (leftmost) bit position, then 
if Q - I. you first shift the previous 
CRC to the right as in figure ll. In 
terms of a high-level language you 
shift the CRC right by one. discarding 
the LSB. and XOR with 8408 By 



testing Q first and then doing the shift 
before the XOR. you avoid the need 
for a 17-bit register. Because of the 
way in which the CRC is stored in the 
variable, the XOR is done with 8408 
rather than 1021. as you might 
expect. This follows directly from the 
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Listing I: CC1TT.C. the CCITT routine for calculating CRCs in C source code. 

BLKCRC(bufptr, crcr.s. counn"* * 1821 M * rev ° rder */ 
unsigned char *bufptr; 
^unsigned int »crcr««, count; 

i /» end BLKCRC */ 

bytecrc(bufptr. crcres) 
unsigned char *bufptr; 
^unsigned int *crcres; 

unsigned int j.ch.Q; 

cT: C ch r > e M.* CrCr9S>>1: 'f* "a XOR % 

/» move next data bit */ 

! /« end bytecrc */ 
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diagram shown in figure 10. Generally 
speaking, if you process data LSB first, 
you can store the CRC with its MSB in 
the MSB position of the integer variable 

Listing 1 contains the implementa- 
tion of a CRC calculation based on the 
CC1TT polynomial. CCITT.C provides 
lots of intermediate results, for clari- 
ty rather than performance. This pro- 
gram implements the CCITT/IBM FCS 
calculation in a standard way with one 
exception: initializing the CRC \Editor s 
note: CCITT.C. XMODEM.C. SDLCASM. 
and XMODEM. ASM /or OS9 are avail- 
able in several formats: see the insert caret after 
page 368. | 

The purpose, however, is to imple- 
ment the de facto XMODEM stan- 
dard. 1 found no published XMODEM 
CRC specification as such. What has 
been published is a C program. 



XMODEM.C. that does the calcula- 
tions (see listing 2), Apart from using 
the CC1TT polynomial, this program 
is not CCITT standard in that data is 
processed MSB first, rather than LSB 
first, and the CRC is initialized to zeros 
rather than to ones. 

Choices in CRC Calculation 
and Transmission 

A CRC calculated by these methods 
has some very desirable error-detec- 
tion properties, but it is not perfect. 
For example, you can add or delete 
any number of zero bits to or from the 
beginning of the block without affect- 
ing the CRC. (The remainder stays 
zero no matter how many zero bits 
start the block.) Furthermore, since 
the CRC is a cyclic code, any error, 
such as a clock slippage, that deletes 



a bit at the beginning of the block and 
inserts the same bit at the end of the 
block (the last bit of the CRC) will not 
affect the CRC. For these reasons the 
CRCs used by IBM in the SDLC pro- 
tocol and CCITT in the X.25 and 
HDLC protocols specify the following: 

1 . All bits of a block are protected 
by the CRC. 

2. The data is sent LSB first. The 
CRC is calculated on bits as they 
are sent. 

3. The CRC is initialized to all ones. 
This allows detection of any 
missed or inserted zero bits at 
the beginning of a block. 
(Missed or inserted ones are still 
detected.) 

4. The one's complement of the 

{[Mlinued) 



Listing 2: XMODEM.C. the XMODEM routine for calculating CRCs in C source code. 
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Sample bit-oriented CRC routine 
Adapted from YMODEM protocol reference 



/* Calculate CRC on a block of data. 

/* Ptr points to block of characters, count gives size of buffer. 
/* This program returns the CRC with the LSB of the CRC in the 
/* high bit of the result integer. 

/* XMODEM deviates from the CCITT standard in that it does not use 
/* the LSB of the data 1st, nor does it initialize the CRC to all 
/* ones as specified by the standard. 

int calcrcfptr, count) 
char *ptr; 
int count; 

i 

unsigned int crc; 
i nt i ; 

crc = 6; 

while ( — count >= 8) \ 
i = (int) *ptr++; 
i - i « 8; 
crc « crc * i ; 



*/ 
*/ 
»/ 
*/ 
*/ 

% 



for (i-0; i<8; ++l) 
if (crc 4 0x8000) i 
crc - (crc « 1); 

crc « crc * 0x1021 ; 



\ /* endif */ 
else 

crc - crc « 1; 

| /* end wh i I e */ 
return (crc 4 0xFFFF); 
\ /* end calcrc */ 



/» note not BxFFFF 

/* convert data char to int 
/* shift char to high byte 
/* add current data to current 
/* remainder modifies only least 
/* sig 8 bits (high byte) of CRC 
/* loop for each bit 
/* test D XOR R0 
/» discard LSB of CRC and 
/* append zero 
/* XOR with low 16 bits 
/* of CCITT polynomial 
/* because CRC is stored LSB 1st 
/« polynomial written MSB 1st 



/* discard LSB 4 append */ 



•/ 

*/ 
*/ 
»/ 

I 



/* 16-bit result for whole block »/ 
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Figure 12: Tne contents of the shift register at the beginning of bytewise SDLC 
calculations. 
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Figure 13: Tne contents of the shift register after the first shift. Note that all 
entries in a column are XORed together. 
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Figure 14: Vie same as figure 13 but with the abbreviation TO for DO XOR RO. 
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Figure 15: The contents of the shift register after three more right shifts. TO 
represents DO XOR RO. T1 represents D1 XOR R1 . and so on. 
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Figure 16: The contents of the shift register after the fifth shift. 
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Figure 17: Tne contents of the shift register after the eighth shift. 



CRC is transmitted rather than 
the CRC itself. This allows detec- 
tion of slippage-type errors. 

5. The CRC is sent LSB first. 

6. The polynomial used is X" + 
X" + X s + 1. 

XMODEM.C implements the CRC as 
follows: 

!. Only data bits are included. The 
header character and two block- 
number bytes are not included. 

2. The data is sent LSB first. The 
CRC is calculated on the data 
MSB first. 

3. The CRC is initialized to zeros. 

4. The CRC is not complemented 
before transmission. 

5. The order of transmission is high 
byte of CRC then low byte. Since 
UARTs transmit LSB first, this is 
equivalent to LSB of high byte 
through MSB of high byte, then 
LSB of low byte through MSB of 
low byte. 

6. The polynomial used is X" + 
X" + X' + 1. 

The CCITT method is easily imple- 
mented in a chip using shift registers 
and XOR gates, while the XMODEM 
method is not easily realized in hard- 
ware. Tb check an incoming data 
block, you have two options. You can 
calculate the CRC on all the protected 
bits only, omitting the CRC bits, and 
compare the calculated value to the 
received value. Or you can calculate 
the CRC on all the protected bits and 
the CRC itself and then compare the 
result to a known constant. If the first 
method is adopted in the CCITT case, 
the one's complement of the calcu- 
lated value must be compared to the 
received CRC. In the XMODEM case, 
the comparison is direct. If the second 
method is adopted, in the XMODEM 
case the known constant is zero, while 
in the CCITT case it is F0B8. (The 
high bit is on the left, i.e., 1.) No one's- 
complementing is required. 

Byte-oriented Software 
Implementations 

The next step is to derive routines to 
calculate the CRC a whole byte at a 
time rather than bit by bit. This ap- 
proach was first proposed by Perez. 
Wismer. and Becker. The motivation 
is that an 8-bit microprocessor is not 
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limited to single-bit XORs but can do 
them 8 bits at a time. The basic ap- 
proach is to work from the hardware 
diagram (figure 9 or 10) and see what 
the CRC register would look like after 
8 bits have been calculated. 

Bytewise SDLC calculations. In 
figure 10 you can see that you have 
the contents of figure 12 in the shift 
register at the beginning of calcula- 
tions. Then you take the LSB of data. 
DO. and XOR it with RO (DO XOR RO). 
After the right shift, the new R15 is DO 
XOR RO, the new R10 is R11 XOR DO 
XOR RO. and the new R3 is R4 XOR 
DO XOR RO. (See figure 13). 

The combinations DO XOR RO, D1 
XOR R1, D2 XOR R2. and so on. oc- 
cur frequently, so let's abbreviate 
them as TO = DO XOR RO, T1 = D1 
XOR RT, and so on. Figure 13 can 
now be rewritten as in figure 14. If you 
proceed in this fashion for three more 
shifts, you get figure 15. 

To do the XOR on D4, you must use 
the content of the LSB of the shift reg- 
ister, which is now R4 XOR DO XOR 
RO, The result is D4 XOR R4 XOR DO 
XOR RO, or in shorter form, T4 XOR 
TO, The result after the fifth shift is 
shown in figure 16: after the eighth 
shift, in figure 17. 

The tedious part is done. Now you 
want to write a program that will pro- 
duce the same result when you're 
working with a byte of data as you 
would get from the shift register after 
eight shifts. The emphasis in this 
routine will be on speed. If you make 
it too general-purpose, a different 
choice of polynomial will lead to a 
completely different program. Also, 
for speed, it makes sense to code the 
routine in assembly language. 

When working with 8-bit microcom- 
puters, it is convenient to define the 
8-bit quantities found in table 2. If you 
study figure 17, you will see that the 
combination 

T7 T6 T5 T4 
XOR T3 T2 T1 TO 

occurs several times. Let's call this 
term U = U7 U6 U5 U4. If you re- 
write figure 17 with these substitu- 
tions, you get figure 18. 

Further optimizations become ap- 
parent as you write the code. It is con- 
venient to implement the SDLC algo- 
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Figure 18: The same as figure 17 but using the abbreviations given in table 2. 
Line 1 is CRCHi moved into CRCLo: line 2 is the high nubble of U and the low 
nibble of T; line 3 is the line 2 byte shifted left by 3 bits: and line 4 is U shifted 
right by 4 bits. 



Table 2: Some convenient abbreviations for various 8-bit quantities used in CRC 
calculations in an 8-bit microcomputer environment. 
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Table 3: Some test cases provided for comparison if you want to write 
your own routine. 

Text SDLC CRC 



T 1B26 14A1 

THE 44BE 7D8D 

THE,QUICK,BROWN 1 FOX,0123456789 DF91 7DC5 
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Figure 1 9: The contents of the shift register at the beginning of bytewise 
XMODEM calculations. 
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Figure 20: The contents of the shift register after the XOR witn the data byte. 



rithm las in the program SDLC.ASM) 
because of the built-in check that cal- 
culating a CRC on a "received" block 
provides: that is. the result should 
always be F0B8. 
The innermost loop of SDLC.ASM 



takes only 86 cycles per byte in- 
cluding seven overhead cycles to 
check for end of buffer. If a 2-MHz 
6809 were dedicated to CRC calcula- 
tions, this would result in a through- 
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CALCULATING CRCS 
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Figure 21: The same as figure 20 6ul vtf the abbreviations T7 for R15 XOR 
D7. T6 for R14 XOR D6. and so on. 
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Figure 22: Ike contents of the s 
performing the left shift. 


$ 1 


■«;$!, 


r after processing the first data bit and 




15 14 13 12 11 10 


9 


8 


7 6 5 4 3 2 1 


T3 T2 T1 TO R7 R6 
T7 T6 T5 T4 


H5 


R4 

T7 


R3 R2 R1 RO T7 T6 T5 T4 
T6 T5 T4 



Figure 23: The contents of the shift register after four shifts. 
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Figure 24: The contents of the shift register after five shifts. 
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Figure 25: The contents of the shift register after eight shifts. 



Table 4: The same test cases as in table 3 but for XMODEM 
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put of about 200.000 bps. A bitwise 
algorithm to do the same job takes on 
average 353 cycles per byte and could 
take as many as 385. So the bytewise 
algorithm is about four times faster. 

The test cases in table 3 are pro- 
vided in case you want to write your 
own routine. The column labeled 
SDUC results when you initialize the 
CRC with ones (this is the result 
before the one's-complementingl. and 
the column labeled CRC results when 
you initialize with zeros. 

XMODEM byte-oriented algo- 
rithms. The XMODEM CRC spec' 
such as it is, is shown in listing 2. If 
you repeat the work done for SDLC 
calculations (in figures 12 through I SI 
using listing 2 as a base, you will ob- 
tain the contents of the diagrams in 
figures 19 through 2 5. The contents 
of the starting CRC are shown in figure 
19. After the XOR with the data byte, 
you get the result shown in figure 20. 
When you abbreviate R15 XOR D7 as 
T7, F/14 XOR D6 as T6. and so on. 
figure 20 becomes figure 21. After 
processing the first data bit and the 
left shift, you get the result shown in 
figure 22: after four shifts, figure 23: 
after five shifts, figure 24: and after 
eight shifts (the final resultl. figure 2 5. 
Next, you define the appropriate vari- 
ables such as CRCHi. CRCLo, T, and 
U. and you look for factors. 

An implementation of the XMODEM 
CRC calculations for the 6809 is given 
in the program XMODEM. ASM. The 
test case results using the XMODEM 
algorithm are shown in table 4. 

Conclusion 

Once you arrive at the basic idea of 
using the XOR function to implement 
modulo 2 division, it is not difficult to 
apply the CRC algorithms, only 
tedious. As with so many other com- 
munications standards, the options 
and the details provide most of the 
difficulties and confusion ■ 

BIBLIOGRAPHY 

Forsberg. Chuck, ed "XMODEM 
YMODEM Protocol Reference Godzilla 
BBS. 15031 621-3746. 

IBM SDLC General Information Manual 
Appendix B. 

McNamara. I E. Tectmwl Asveas of Data Com- 
munications Bedford. MA: Digital Equip- 
ment Press 1982. 



- 



SYSTAT 



MORE STAT] 




How has SYSTAT become one of the largest 
statistical software companies in only 3 years? 



LESS BULK: 

■ 

Will run Irom floppy disks 
■ 

Needs less memory (256K on IBM or 
compatible machines, 64K on CP/M. 400K 
on VAX, or 512K on Apple Macintosh) 
■ 

Has fewer than Vs of the commands of 
other manufacturers' packages 

MORE STATISTICS: 
■ 

Full screen spreadsheet data editor 
■ 

Missing data, arrays, character variables 
■ 

Unlimited cases 
■ 

Process rectangular, hierarchical, triangular 

files and variable length records 

■ 

Relational database management and 

tile concatenation 

• 

Character, numeric, and nested sons 
Unlimited numeric and character 



Interactive or batch 
■ 

Read and write text and external files 
■ 

Subgroup processing in statistical modules 

with SELECT and BY 

■ 

Value labels 
■ 

RECODE statements for quick multiple codes 

■ 

Scatterplots. contours, histograms, stem- 
and-leaf. boxplots, bar chans, quantiie, 
probability plots 
■ 

Basic statistics, frequencies, T-tests 
■ 

Multi-way crosstabs with log-linear modeling, 
association coefficients. PRE statistics, 
asymptotic standard errors 
■ 

Pairwiselistwise missing value correlation. 
SSCP, covanance, Spearman. Gamma, 
Kendall Tau. Euclidean distances 
i 

Linear, polynomial, multiple, stepwise. 

weighted regression 

■ 

Extended regression diagnostics 
■ 

Multivariate general linear model 



Multi-way ANOVA. ANOCOVA. MANOVA, 
repeated measures, unbalanced designs, 
post-hoc tests 
■ 

Principal components with rotations 
and s 



Multidimensional scaling 
■ 

Multiple and canonical discriminant analysis. 

Bayesian classification 

■ 

Canonical correlation 
■ 

Cluster analysis (hierarchical, single, 
average, complete median, centroid 
linkage, k-means, cases, variables) 
■ 

Nonparametric statistics (sign, Runs, 
Wilcoxon, Kruskal-Wallis, Friedman two-way 
ANOVA, Mann-Whitney U. Kolmogorov- 
ltof 



Time Series (smoothing, seasonal and 
nonseasonal ARIMA. ACF, PACF, Cross- 
correlation function, transformations. 



Nonlinear estimation (non-linear regression, 
least absolute values regression, logit, probit. 
maximum likelihood estimation, iteratively 
reweighted least squares) 
■ 

Supplements include: Logit (multinomial 
logit). Probit (probit analysis). Testal (classi- 
cal and Rasch model test analysis), Report 
(formatter tor SYSTAT files), Transter (dBase, 
Lotus, GAUSS, STATA, SPSS direct file 
transfer) 

SYSTAT operates on the following machines: 
IBM-PC XT AT.- Apple II,' - Apple Macin- 
tosh,'" Kaypro.'" HP 1 50. " HP 9000. " DEC 
Rainbow.'" VAX.'" Alpha Micro.'" MS-DOS- 
CP M'" and UNIX.'" 

Single copy price: 

S595 USA and Canada ' $695 Foreign 
Site licenses and quantity prices available 
Call or write for additional information 
SYSTAT, Inc. 
2902 Central Street 
Evanston. IL 60201 
312 864.5670 
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Read more about the extraordinary success 
ot SYSTAT in So You ve Go! a Great lOea 
by Steve Filler (Addison Wesley. 1986). 
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